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Abstract. We carefully reanalyze the ROSAT PSPC X- 
ray spectro-photometric observations of HCG 16 (Arp 
318), and compare them to optical and radio data. Its X- 
ray morphology resembles its morphology at 20 cm, seen 
by the NVSS. In particular, we detect diffuse emission in 
eight regions filling half of the 200/i 50 1 kpc (8.'7) radius 
circle around the optical center of the group: one region 
encompassing galaxies a & b, two regions surrounding the 
group galaxies c & d, a clumpy region roug hly 140 kpc 
from the group galaxies, which may be gas ejected from 
one of the galaxies, plus regions respectively associated 
with a background radio-source, a probable background 
radio-source, a foreground star and a background group 
or cluster. The bolometric X-ray luminosity of the dif- 
fuse emission, excluding the regions associated with ra- 
dio galaxies, is L^°' = 2.3 x 10 41 /i^ 2 ergs _1 , i.e., half of 
the luminosity found by Ponman et al. (1996). The region 
that is offset from the galaxies contributes half of the dif- 
fuse X-ray luminosity of the group. The diffuse emission 
is cool (T < 0.55 keV with 90% confidence with a best fit 
T = 0.27 keV). At these low temperatures, the correction 
for photoelectric absorption in the estimate of bolometric 
luminosity is a factor 3.5 and varies rapidly with temper- 
ature, hence an uncertain bolometric luminosity. 

The clumpy distribution of hot diffuse gas in HCG 16 
is illustrated by the low mean X-ray surface brightness 
and hot gas density of the regions of undetected emission 
within 8' (at most 1 /4 and 1 /6 of those of the detected 
gas, assuming both have same temperature, metallicity 
and dumpiness). The irregular X-ray morphology of the 
diffuse emission rules out a (nearly) virialized nature for 
HCG 16, unless intergalactic gas had sufficiently high spe- 
cific entropy to be unable to collapse with the group. In 
any event, the clumpy gas distribution, and high luminos- 
ity given the low temperature suggest that most of the 
diffuse gas originates from galaxies, either through tidal 
stripping or through galactic winds driven by supernova 
remnants. Therefore, no spiral-only HCGs are known with 
regular diffuse emission tracing a gravitational potential. 
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Our results highlight the need for a careful 2D spa- 
tial analysis and multi-wavelength study of the diffuse X- 
ray emission from groups, suggesting that other compact 
groups could be significantly contaminated by superim- 
posed X-ray sources. 
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1. Introduction 

The extreme apparent density of compact groups of galax- 
ies makes them ideal sites for the study of multiple in- 
teractions of galaxies, and more generally, another dense 
environment to be compared with rich clusters. Indeed, 
compact groups (of typically 4-5 bright galaxies) such as 
those (hereafter, HCGs) cataloged by Hickson (1982) ap- 
pear as compact in projection onto the plane of the sky 
as the cores of rich clusters, and are moreover selected to 
be isolated. 

The discovery of numerous signs of galaxy-galaxy in- 
teraction within HCGs (see Hickson 1997 for a review) 
suggest that most HCGs are indeed dense in 3D. However, 
the very short crossing times derived from galaxy spec- 
troscopy (Hickson et al. 1992) suggest rapid galaxy merg- 
ing and coalescence into a single giant elliptical galaxy 
(Carnevali et al. 1981; Barnes 1985; Mamon 1987; Barnes 
1989; Bode et al. 1994; Governato et al. 1996) if indeed 
they correspond to the low-mass end of clusters of galax- 
ies, forming at high redshift, unless the group is constantly 
replenished through accretion of surrounding galaxies (Di- 
aferio et al. 1994; Governato et al. 1996). In fact, a variety 
of arguments have been put forth suggesting that compact 
groups are mostly chance alignments of galaxies along the 
line of sight within larger systems: loose groups (Rose 
1977 for chain-like groups; Mamon 1986 for the major- 
ity of compact groups), clusters (Walke & Mamon 1989) 
or cosmological filaments (Hernquist, Katz, & Weinberg 
1995). Thanks to gravity, these chance alignments tend 
themselves to be binary-rich (Mamon 1992), and it is very 
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difficult to tell whether the interaction seen in HCGs are 
caused by a system of 4 or more bright galaxies or simply 
by binaries, well-separated along the line-of-sight. 

Recently, there has been much hope that the debate 
on the nature of HCGs could be resolved by X-ray obser- 
vations. The ROSAT, ASCA and Beppo-SAX X-ray satel- 
lites are sensitive enough in the soft X-ray band to be able 
to detect the diffuse intergalactic plasma within (nearly) 
virialized galaxy systems with potential wells with depth 
corresponding to a ID velocity dispersion of ~ 250 km s _1 . 
Most sensitive of these is the Position Sensitive Proportional 
Counter (PSPC) of ROSAT, which discovered diffuse inter- 
galactic emission from large numbers of compact groups 
(Ebcling, Voges, & Bohringer 1994; Pildis, Bregman, & 
Evrard 1995; Saracco & Ciliegi 1995; Mulchaey et al. 1996; 
Ponman et al. 1996, hereafter PBEB). 

But there has been a debate on whether the fraction 
of compact groups detected by ROSAT is 40% (Mulchaey 
et al. 1996) or 75% (PBEB, taking into account the selec- 
tion effects on distance) . There is also controversy on the 
luminosity-temperature relation: Mulchaey & Zabludoff 
(1998) derive a relation that is consistent with the extrap- 
olation from rich clusters, while PBEB find a luminosity- 
temperature relation for compact groups with a much 
higher slope, with the hottest groups lying on the clus- 
ter extrapolation. 

Moreover, the morphology of the diffuse X-ray emis- 
sion of compact groups is very diverse, as attested by 
the PSPC maps provided by Pildis et al. (1995), Saracco 
& Ciliegi (1995), Ponman & Bertram (1993) (only for 
HCG62) and PBEB (only for HCG 16). In some cases, 
such as HCG 62 (Ponman & Bertram), diffuse intergalac- 
tic emission extends well beyond the group, not centered 
on any galaxy of the group, with regular circular isophotes, 
just as is to be expected in a well relaxed galaxy system. 
But in other cases, the emission is only attached to in- 
dividual galaxies, see e.g. HCG 44 (PBEB). And there 
are intermediate cases, such as HCG 16, where the dif- 
fuse emission does not appear as extended as the galaxy 
system nor as regular as in HCG 62. 

In fact, HCG 16 (also known as Arp 318) is an unusual 
galaxy system. First, because the six brightest galaxies of 
the group (Hickson's original 4 plus two more outside the 
group isolation annulus, see de Carvalho et al. 1994) are 
starburst, LINERs or AGNs (Ribeiro et al. 1996, here- 
after RdC3Z). Moreover, the X-ray properties of HCG 16 
are controversial and possibly extreme. It was first de- 
tected with the EINSTEIN satellite (Bahcall et al. 1984), 
which did not have the angular resolution to resolve the 
emission between the group galaxies and an intergalactic 
medium. However, Ponman et al.'s analysis of ROSAT 
PSPC observations made it the coldest detected group 
(T = 0.30 ± 0.05 keV), and there are no other spiral-only 
compact groups with diffuse X-ray emission (Mulchaey 
1999; see also Ponman et al.). HCG 16 is thus an abnormal 
group given the very strong correlations between X-ray lu- 



minosity and group spiral fraction found by Pildis et al. 
(1995) and Mulchaey et al. (1996). Moreover, whereas dif- 
fuse X-rays were clearly detected by PBEB, Saracco & 
Ciliegi (1995) failed to detect such diffuse emission at an 
upper limit 16 times lower, 1 whereas only a factor 2.3 
(which we find by simulating a MEKAL plasma with tem- 
perature, abundance and absorbing column as quoted by 
PBEB) is attributable to the wider ( "bolometric" ) energy 
range in which PBEB compute their luminosities. Given 
the low temperature that PBEB derive for HCG 16, their 
derived X-ray luminosity places it two orders of magni- 
tude above their compact group luminosity-temperature 
relation and roughly a factor of two above the extrapola- 
tion of the cluster trend. It thus seems difficult to reconcile 
HCG 16 with a low temperature extrapolation of regular 
X-ray emitting compact groups. Indeed, PBEB note that 
HCG 16 is "probably not fully virialized" . 

In this article, we present a detailed analysis of the 
ROSAT PSPC observations of HCG 16. The data reduc- 
tion is presented in Sect. 2, our spatial analysis in Sect. 3, 
and our spectral analysis in Sect. 4. In Sect. 5, we com- 
pare our results with previous X-ray analyses of HCG 16, 
perform a mass budget of the group and ask if it is viri- 
alized. In a following paper (Mamon & Dos Santos 1999, 
hereafter Paper II), we discuss at length the cosmological 
and dynamical constraints on the nature of HCG 16. 

2. Observations and data reduction 

2.1. Observations and preliminary reduction 

HCG 16 was observed in January 1992 with the PSPC 
(in its low-gain state), on board the ROSAT satellite, 
for a total observing time of 14 634 s. We obtained the 
data from the ROSAT archives. Snowden et al. (1994)'s 
PSPC Extended Source Cookbook software was used to per- 
form the first-pass data reduction, i.e., rejection of high- 
background times, energy-dependent (in 7 bands) back- 
ground subtraction, exposure and vignetting corrections. 
We adopted a conservative value of 170ctss _1 for the 
maximum Average Master Veto rate allowed (see Snow- 
den et al. 1994). Even with this low threshold, only 6% 
of the total observing time was rejected, leaving an ef- 
fective observational time of 13 748 s. We then carefully 
examined the light curves of the total counts in the en- 
tire image per energy band as defined by Snowden et al. 
(1994), and checked that no short time scale glitches were 
present. Point sources were detected using a sliding box 
algorithm, with the improvement that the box is a cir- 
cle with a radius varying with off-axis angle, to model 
the varying point spread function. Each point source de- 
tected at a level exceeding 3cr was removed, i.e. a circle 

1 Given the fluxes measured by Saracco & Ciliegi (1995) for 
HCG 16 and their adopted value for H , their quoted upper 
limit for their luminosities were underestimated by a factor 2 
for all undetected groups in their Table 4 except HCG 3. 
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centered on the source, of radius 1.5 times the 90 % encir- 
cled energy radius, was excised. Unless otherwise stated, 
we limited our analysis to the [0.2, 2.0] keV energy band, 
because at lower energies, the background is too high to 
be adequately subtracted from the data, while at high en- 
ergies the same occurs because our sources happen to be 
relatively cold, and moreover the calibration of the PSPC 
is uncertain. This preliminary reduction produced an im- 
age with 512 x 512 pixels, 15" wide (roughly the FWHM 
of the PSPC's PSF at IkeV). 

2.2. Background estimation 

The definition of the background region, whose counts are 
subtracted to each pixel within some region, is of crucial 
importance for the spatial detection of sources with very 
low signal-to-noise ratio (hereafter S/N), as well as for 
spectral analyses. 

The shadowing by the supporting ring of the PSPC 
(situated at ~ 20' from the center of the field) is visible 
in the images, even after the vignetting correction. Thus, 
we measure the background well outside of the ring. In 
practice, we choose three annular regions to measure the 
background, each centered on HCG 16, with inner and 
outer radii of 30' - 48' (BG 1 ) , 26' - 40' (BG2) , and 26' - 34' 
(BG3). The radial structures of the PSPC supporting ring 
were removed in each case. Table 1 shows the background 
counts in the three regions within the 0.2 — 2.0 keV energy 
range. We note that BG3 has a slightly lower value than 
the other two. 

Table 1. Background from different annular regions 



Region 


Radii 


Pixels 


Counts/pixel 


BG1 


30' - 48' 


25930 


0.593 ± 0.005 


BG2 


26' - 40' 


17630 


0.577 ± 0.006 


BG3 


26' - 34' 


9832 


0.565 ± 0.008 



To decide which is the best background, we measure 
the net counts in two regions within the inner ring sup- 
port, using each of the three background regions for mea- 
suring the background. The latter have been vignetting 
corrected and normalized to the number of pixels of each 
region. Our two test regions are the 8' radius circle cen- 
tered on the group optical center and the annulus sur- 
rounding this circle, with inner and outer radii of 8' and 
17'. Both arc in the inner 20' of the field of view, thus 
avoiding problems with the supporting ring of the PSPC 
and with somewhat uncertain vignetting correction out- 
side of this ring. Table 2 gives the background subtracted 
data and errors for both regions and for each background. 

While we find positive net counts within the inner 8' of 
HCG 16, whatever three of the background regions is used 
to estimate the background, only BG3 is compatible with 
the counts within the 8' — 17' annulus, yielding near zero 



Table 2. Tests of background subtraction 

Region r < 8' 8' < r < 17' 

Net counts Significance Net counts Significance 

BG1 222 ± 48 4.6 a -250 ± 78 -3.2 a 

BG2 271 ±49 5.5 a -121 ± 82 -1.5 a 

BG3 308 ±51 6.0 a -26 ± 91 -0.3 a 



net counts. The net counts with the other two background 
regions are difficult to understand, unless there happens 
to be X-ray absorption by Galactic or intergalactic neutral 
hydrogen merely in this 8' — 17' annulus. This appears to 
be ruled out by HI observations with the VLA (Williams 
1998; Verdes-Montenegro 1999) and by the spectral anal- 
ysis (Sect. 4). We thus infer that BG1 and BG2 are con- 
taminated by sources or suffer from uncertain (large) vi- 
gnetting corrections. In fact, Pildis ct al. (1995) had also 
encountered a PSPC field with a background that rose 
with distance to the field center. We therefore use BG3 
(26' < r < 34') to measure the background, which then 
amounts to B = 6.6 ± 0.1 x 10 -4 s _1 arcmin~ 2 (in the 
[0.2 — 2.0] keV energy range). 

3. Spatial analysis 

3.1. Preliminary spatial analysis 

Ponman et al. (1996) claimed a 3cr detection of diffuse 
emission within a radius ~ 8.'5 surrounding HCG 16 (cor- 
responding to ~ 195 Ii^q kpc, given HCG 16's redshift of 
0.0132 2 ). The net counts shown in Table 2 confirm PBEB's 
global detection of diffuse emission in HCG 16. We detect 
diffuse emission at 6 a within 8' and fail to detect sig- 
nificant counts between 8' and 17'. Hence, to first order, 
the extent of the diffuse emission is very roughly 8', or 
190/i 50 1 kpc, similar to what is inferred from PBEB's sur- 
face brightness profile of HCG 16. 

But we wish to go further: what is the spatial distribu- 
tion of this excess of photons within 8' of the group center? 
Here we are confronted with the low number of detected 
photons (~ 300). We applied three different methods to 
analyze the photon distribution within 8'. First, we count 
the photons in a grid encompassing the 8' radius circle. 
Then, in order to obtain the best S/N ratio, we adaptively 
smooth the image and detect diffuse emission. Finally, we 
apply a wavelet-based method to detect structures at all 
scales and verify that emission is present on more than 
one scale as a confirmation of its diffuse nature. 

3.2. Analysis on a grid 

The simplest way to spatially analyze HCG 16 is to define 
a grid encompassing the whole group. We divide the field 

2 Distances throughout this paper are derived assuming a 
Hubble constant Ho — 50 kms -1 Mpc -1 , with hm = 1, 
unless explicitly given. 
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Fig. 1. Contour map of smoothed (by a Gaussian of FWHM 
= 45") ROSAT-PSPC X-ray emission from HCG 16 super- 
imposed on an optical image of the group provided by the 
Digital Sky Survey (DSS). Coordinates are for epoch J2000. 
This figure and the following ones do not mask the de- 
tected point sources (which are masked in the subsequent X- 
ray photometry and spectral analyses). Contours are drawn 
at 1 a, 2 a, 3 a, 4 a, 5 a, 10 a, 15 a, 20 a, 30 a, and 50 a above the 
background level, where a is the standard deviation of the 
smoothed background. The scale (upper left) is for the plane at 
the distance of HCG 16. The four original galaxies discovered 
by Hickson (1982) are marked 'a', 'b', 'c' and 'd', together with 
galaxy HCG 16-3 (de Carvalho et al. 1994), marked '3', which 
is at the same distance as the group. The dotted circle centered 
on the group optical center is 8.' 7 wide, i.e., 200/i^ 1 kpc at the 
distance of HCG 16. The dashed square grid used for the pre- 
liminary spatial analysis has cells 16 pixels wide (92 /i^kpc). 
Cells are counted from left to right then from top to bottom. 



of view of the group in sixteen 16 x 16 pixel squares, (i.e., 
92 ft-g^kpc side at the distance of HCG 16). The grid over- 
laid on a smoothed X-ray image of HCG 16 is displayed in 
Fig. 1. After removing the point sources (including group 
galaxies), we counted the background-subtracted number 
of photons in each square of the grid. The results, as well 
as the statistical significance of the detections, are given 
in Table 3. 

The analysis of Table 3 shows that 7 regions over the 
16 selected are detected at a 2 a level above the back- 
ground. We count regions from East (left) to West (right) 
and North (top) to South (bottom). Fig. 1 shows that 
the regions of excess counts are concentrated East, West 
and South of the four bright galaxies of the group: regions 
9, 13 and 14 in the East contain 114 net counts, while 
regions 8, 11, and 12, contain 149 net counts. Moreover, 



Table 3. Net counts within a square grid 



Region 


Pixels 


Net counts 


Significance (it) 


1 


106 


-4 ± 9 


-0.5 


2 


252 


-9 ± 13 


-0.7 


3 


204 


-3 ± 12 


-0.2 


4 


217 


17 ± 13 


1.3 


5 


256 


-8 ± 13 


-0.6 


6 


198 


23 ± 13 


1.8 


7 


167 


31 ± 13 


2.5 


8 


256 


57 ± 16 


3.7 


9 


200 


50 ± 14 


3.6 


10 


203 


15 ±13 


1.2 


11 


194 


34 ± 13 


2.5 


12 


256 


58 ± 16 


3.8 


13 


200 


36 ± 14 


2.7 


14 


256 


28 ±15 


2.0 


15 


196 


18 ±13 


1.4 


16 


201 


-18 ±11 


-1.6 



The regions, shown in Fig. 1, are numbered from East to West, 
looping North to South. Cells with a statistical significance of 
detection greater than 2 a are shown in bold. Errors are 1 a 
and on the total (background + net) counts. 



half of the cells in the grid (if we don't take into account 
cell 6, detected at a 1.8 a level) are compatible with no 
excess emission over the background, and most of these 
are North of the four galaxies (regions 1, 2, 3, 4, 5, 10, 
15 and 16). Hence, the diffuse emission is located in two 
distinct regions as well as around the galaxies. This could 
be seen in PBEB's Fig. 2, but they gave no quantitative 
analysis of the spatial distribution of photons, except for 
their radial surface brightness profile. 

Is it possible to detect more precisely these excess 
counts regions, without degrading the S/N ratio? Indeed, 
the grid we used was arbitrarily set on the HCG 16 image, 
and some cells, especially cells 10, 11, 13 and 14, over- 
lap two regions where the photon densities differ. Using a 
smaller grid would not help because we are at the limit of 
sensitivity. Since we are principally interested in mapping 
the diffuse intergalactic gas, we can smooth the image, 
looking for large-scale features. But we need to be care- 
ful with the level of smoothing: indeed, regions with low 
count numbers must be smoothed on larger scales than 
the bright regions (e.g. , the bright galaxies), so as to 
keep a good statistical significance of the regions of diffuse 
emission. Consequently, we choose an adaptive filter algo- 
rithm, which automatically adapts the smoothing length 
to the local density of photons. PBEB already used adap- 
tive smoothing to detect diffuse gas in HCG 16 (see their 
Fig. 2), but did not make a quantitative use of the infor- 
mation obtained with this technique. 
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3. 3. Adaptive filtering of the image 
3.3.1. Spatial resolution 

In a classical (top-hat) smoothing filter, the filtering radius 
is fixed, the smoothed intensity is the mean (unweighted) 
counts within this radius, and the total counts within this 
radius varies across the image. Hence, every pixel will have 
a S/N directly proportional to the local photon density. 
The pixels in low surface brightness regions will then have 
poor statistics. 

In an adaptive filter, the filter size is adapted to contain 
a fixed number of counts. Let Co be this fixed number of 
(background + net) counts per smoothing region of size P 
pixels. The intensity of a given pixel will typically be I = 
C /P, hence S/N = I /SI = P/5P. Now, if the cumulative 
counts rise with radius 1 a faster than on average, they 
will reach Co at P — SP, where the typical counts are 
Ci = I(P — SP). Since the count curve of growth, C(P), 

1/2 

is a Poisson process, one has Co = C\ +C X . Then, in the 
limit SP < P, one obtains SP = (P//) 1/2 , hence S/N = 
C 1/2 (we check that SP/P = C~ 1/2 < 1). Therefore, 
every pixel will have the same S/N. 

Refining the grid analysis, we define polygonal regions 
with much greater S/N than in the grid cells. Contrary to 
the grid, there will be no regions half overdense and half- 
underdense in photons. We produce adaptively smoothed 
images using ADAPT in the PSPC Extended Source Cook- 
book (Snowden et al. 1994), with C = 25, 50 and 100. 
Before smoothing, we masked point sources detected at 
the 3 a level, so as to better highlight the diffuse emission. 
The adaptively smoothed images are shown in Fig. 2. For 
esthetical reasons, the point sources detected by DETECT 
are not cut in these images. 

Our smoothing resolution is easily computed. Calling 
B the mean background, since I > B, our smoothing ra- 
dius, equal to (P/tt) 1 / 2 , must be smaller than 

" 2 • (i) 

With a mean background of B — 0.565 counts/pixel (see 
Sect. 2.2), we obtain the smoothing resolutions listed in 
Table 4. 



Table 4. Adaptive smoothing parameters 



Co 
(cts) 








S/N 


(pixels) 


(arcmin) 


(ftgokpc) 


25 


3.8 


0.9 


22 


5.5 


50 


5.3 


1.3 


31 


7.6 


100 


7.5 


1.9 


43 


10.5 



The smoothing radii are obtained through Eq. (1) with a mean 
background of B = 0.565 counts/pixel (see Sect. 2.2). 

The grid we used in Sect. 3.2 was made of 16 x 16 pixel 
cells. Table 4 shows that the lowest resolution varies with 



Co from 1/4 square cell (i? sm ~ 4 pixels) for Co = 25 to 
1 square cell (i? sm ~ 8 pixels) for C = 100. This shows 
that our adopted values of C are well suited for the size 
of the group, and will not smear out intermediate scale 
structures by smoothing on too large a scale. 

3.3.2. Results 

Consider first the image smoothed with C = 100 
(Fig. 2a), which has the largest features of the three im- 
ages. Within the single contour encompassing the group 
galaxies are two overdensities, SE and SW of the group. 
They correspond respectively to part or all of cells 9, 13, 
14 for the SE overdensity and cells 8, 11, 12 and 15 for 
the SW one, in the grid of Sect. 3.2. These cell numbers 
correspond to the cells detected at a level higher than 
2cr (see Table 3). We divide this closed region in 5 dif- 
ferent polygonal regions (C1-C5) filling the contours of 
the X-ray image adaptively smoothed with Co = 100, as 
shown in Fig. 2. The exact cut between different regions 
is arbitrary, but we tried to be consistent with what we 
know about HCG 16: CI is the South-Eastern overdensity, 
C2 and C3 are the regions around galaxies c & d and a 
& b of the group (but with these galaxies cut out), C4 
is the South- Western overdensity, and C5 is a somewhat 
disjointed region, to the West of C4. The source to the 
South of C4 is detected as a point-like X-ray source, and 
also appears as a double point-like radio-source (Fig. 2d) . 
Finally, C6 is the region within the 8' radius circle minus 
the closed region (C1+C2+C3+C4+C5), i.e the part of 
this circle (53% of it) with no obvious excess counts over 
the background. 

Table 5 shows the net counts in each region. Regions 
CI and C4 are detected at ~ 5cr, much higher than re- 
gions C2 and C3 surrounding the bright galaxies (detected 
at ~ 3 a). This shows that the adaptive smoothing has tar- 
geted regions of high S/N ratio more efficiently than the 
spatial detection on the square grid. But the most impor- 
tant feature, which fully justifies the use of adaptive filter- 
ing, is that the emission of region C6 is fully compatible 
with the background. Hence, half of the region within 8' of 
the group center contains negligible emission from diffuse 
hot gas. In Sect. 4.1 below, we confirm the different nature 
of the two halves of the 8' circle. 

3.3.3. Detailed spatial analysis 

Nothing has yet been said about the nature of the X- 
ray photon overdensities detected. We have intentionally 
called them with the vague denomination "regions" . Are 
these photons emitted by diffuse gas linked to the group? 
If so, is this gas primordial, or mainly ejected by the galax- 
ies? Or does the diffuse emission originate from foreground 
or background sources, not necessarily linked with the 
group? Indeed, looking at PBEB's Fig. 2, the position an- 
gle of the northern overdensity seems to be well-correlated 



6 



S. Dos Santos & G. A. Mamon: Clumpy diffuse X-ray emission from HCG 16 



m 



Si 



1|Hh 



02 10 m 30 3 15 9 lO^OO 3 45 5 30 9 15 5 09 m 00 5 45 9 

R.A. 



CD - 

a - 



:;:; Kra^M^ ;:;:;:;:;;:; ^ ;;:;;;: 



■fill 



-.1 '. ■» ■ 



Itt C1C ! 



D 



Si 



02 1Q m 3Q E 15 s 10 m 00 a 45 a 30 s 1 5 a Q9 m Q0 a 45 s 

R.A. 



o 

CD - 

a - 



NliiSEIIiSSIj!: 
I 2= 1' 



jPI ^5 



(i-V 



Ipiii 
iiiiiiii! 



.!■ . cm 

* cm .ni 



02 1Q m 3Q £ 15 s 10 m 00 a 45 a 30 s 1 5 a 09 m 00 a 45 s 

R.A. 



o 



-T-t , 



111 



Biiiii 



IPfP 



XR3 



,1 



ll!l!I|pI:IlIiI!IMp!l!l!!IlIII 
liIlBlSIiJillBil| 



t^ 3 



NVSS 



02 10 m 30 s 15 s 10 m 00 a 45 a 30 s 1 5 a 09 m 00 a 45 £ 

R.A. 



Fig. 2. Contour maps of the adaptively smoothed X-ray emission (Co = 100 (a), 50 (b), and 25 (c), see Table 4 and Eq. 1) 
superimposed on an optical DSS image. The five polygonal regions (C1-C5) dividing the emission region in HCG 16 are also 
shown, as well as the different components of region CI (CIA, C1B, C1C and C1D, see text). Lower right plot (d) shows a 
contour map of a 20 cm radio image from the NVSS (see Condon et al. 1998) . 



with the alignment of three optical sources. The spatial 
and spectral capabilities of ROSAT are certainly insuffi- 
cient to answer these questions, but the careful analysis of 
the images and spectra can provide some useful insight. 
Observations of HCG 16 at other wavelengths can also 
help to specify the dynamical state of the gas. In particu- 
lar, deep optical images and radio surveys can precise the 
interactions between the diffuse gas, the galaxies and the 
radio sources. 



We can now estimate the number of independent 
smoothing regions within our regions of excess counts. If 
our region has area A pixels and E excess counts, the 
number of independent smoothing regions within it is 



TV = 



A AB + E 



Co 



(2) 



Using Table 5 and Eq. (2), we find that, for C = 100 
(our worst spatial resolution), our regions CI, C2, C3, 
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Table 5. Net counts in regions denned after adaptive filtering 



Region 


Pixels 


Net counts 


Significance (a) 


CI 


447 


96 ± 19 


5.1 


C2 


261 


44 ± 14 


3.2 


C3 


206 


45 ± 13 


3.5 


C4 


306 


83 ± 16 


5.2 


C5 


175 


28 ± 11 


2.5 


C6 


1697 


10 ±31 


0.32 



The counts are computed from the unsmoothed images, within 
the 0.2 — 2.0 keV energy interval and are rounded to the nearest 
integer. Errors are 1 a on the total (background + net) counts. 



C4, and C5 consist of 3.6, 2.0, 1.7, 2.7, and 1.3 inde- 
pendent smoothing regions, respectively. Hence, the cor- 
relation among neighboring pixels introduced by adaptive 
smoothing may connect several local maxima in the X-ray 
surface brightness map, which are close in the plane of the 
sky, but not necessarily linked, neither among them nor 
with the group. 

This is illustrated by comparing the adaptively 
smoothed images with different parameters Cq. Compar- 
ing Figs. 2a and 2c, (with Co = 100 and 25 respectively), 
we see that region CI is composed of an elongated struc- 
ture, almost perpendicular to the group (marked CIA), 
together with three point sources (marked C1B, C1C and 
C1D). This substructure remains in remains at Co = 50, 
but is smoothed out in Fig. 2a into the entire CI region. 
All these structures can be seen in Fig. 1, although with 
a worse S/N. On the contrary, region C4 seems to be ex- 
tended, even with the lightest smoothing. 

3.3.4. Spatial positioning of the images 

We have obtained an optical image of HCG 16 from the 
Digitized Sky Survey (DSS), as well as a 20cm contin- 
uum radio image from the NRAO VLA Sky Survey (NVSS), 
which we superpose on a ROSAT PSPC image. The 
FWHMs are ~ 2" (DSS), 30" (PSPC) and 45" (NVSS). 
The radio and X-ray images both have 15" pixels, while 
the optical image has 1.8" pixels. Fig. 2d shows the NVSS 
radio contours overlaid on a DSS optical greyscale image. 
Both optical and radio images were centered at the center 
of the X-ray image. To ensure a correct superposition, we 
aligned the three images using the bright sources detected 
in the three wavebands. Comparison of Figs. 2b, c and 2d 
clearly shows that the X-ray and 20 cm morphologies are 
remarkably similar. 

In their spectroscopic survey of HCG 16, Ribeiro et al. 
(1996) found that six galaxies among the seven belong- 
ing to the dynamical group were emission-line galaxies, 
AGN, LINERs or starburst galaxies. These galaxies are 
good candidates for detection in all three wavebands. In- 
deed, we find four sources detected in the three images: 
they are all HCG 16 galaxies and are marked a, b, c and 



d in Fig. 2d. Galaxies HCG 16a and HCG 16b seem to 
share a common X-ray and radio halo, thus enhancing 
the probability of tidal interaction between these galax- 
ies suggested by optical tails seen in galaxy HCG 16a. 
The differences between the optical, radio and X-ray posi- 
tions are less than 5", much less than the radio and X-ray 
PSFs. Moreover, we find two objects exactly coincident in 
radio and X-rays, ensuring good correspondence between 
X-rays and radio images independently. Both were first 
detected in the X-rays by Saracco & Ciliegi (1995), who 
called them XR1 and XR2, with no optical counterparts. 
We found two radio sources within a radius of less than 5" 
around them, and XR1 has the same bimodal structure in 
X-rays and in radio. 

The probability of finding an NVSS source within a 
given radius around an arbitrary position is less than 10 -3 
when this radius is ~ 5" (see Condon et al. 1998). This 
means that the coincidences between six sources in the 
whole image cannot have occurred by chance, and this en- 
sures that the superposition of the three images is quite 
perfect, given the positioning uncertainties. We are now 
able to compare small-scale structure in the three wave- 
bands and attempt to elucidate the nature of the three 
regions of excess X-ray counts. 

3.4- Finalizing the regions of diffuse emission of HCG 16 

The problem is now to separate local maxima due to in- 
terloping X-ray sources (point or extended sources) from 
those due to the presence of diffuse gas. The regions C2 
and C3, surrounding the four group galaxies, are certainly 
related to the group. We must then analyze the small scale 
X-ray structure of regions CI, C4 and C5 and compare it 
to optical and radio images. For this, we use a wavelet 
transform of the image to search for structures simultane- 
ously at all scales. 

3.4.1. Wavelet transform of the X-ray image 

Since its invention in the early 80s, the wavelet trans- 
form (hereafter, WT) has proven its capabilities in nu- 
merous astronomical applications, such as the detection of 
the large scale structure {e.g., Slezak et al. 1993), galaxy 
detection and counts (Slezak et al. 1990), and structure 
detection in low-intensity X-ray images (Starck & Pierre 
1998). We used the TRANSWAVE a trous (Shensa 1992; 
Starck & Murtagh 1994) wavelet package kindly provided 
by E. Slezak. In the a trous implementation of the dis- 
crete WT, an N x N image is transformed into i wavelet 
planes (hereafter WPs) of N x N pixels, each being the 
difference between two consecutive wavelet smoothings at 
scales i and i + 1 (with 2* and 2 l+1 pixels respectively). 
The pixel values in these planes are called the wavelet co- 
efficients at scale i. The main advantage of this algorithm 
is that each WP has the same number of pixels, and thus, 
the reconstruction of the image (for example, after thresh- 
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olding in the wavelet planes) is a straightforward process 
of addition. However, there is redundancy in the full set of 
wavelet coefficients. The main difficulty with the a trous 
wavelet filtering is the estimation of the statistical signif- 
icance of the pixels in each WP, which should not follow 
Poisson nor gaussian statistics. 

Fig. 3 shows the wavelet contours of WPs 2, 3, 4 and 
5, respectively. We do not take WP1 into account, since it 
is highly contaminated by small-scale noise in the original 
image, nor WP6 and higher, which smooth the group on 
too large a scale to prove useful. Structures seen only at 
the 1 a level are should mostly be artefacts. 

There are two ways we can ensure that a detected 
source is real: First, if a source is detected in more than 
one WP, its probability of being true is enhanced. In- 
deed, random noise is not correlated between WPs and 
cannot produce large spatial overdensities. Moreover, it 
cannot produce pixels high enough to leave some power 
on larger scales. Second, the probability of a false X-ray 
detection being randomly superimposed on an optical or 
radio source is very low. Consequently, we will search for 
small-scale X-ray sources in regions CI, C4 and C5, de- 
tected in several WPs and roughly superimposed on opti- 
cal and/or radio sources. 

3.4.2. Nature of region CI 

The multiple nature of region CI (Sect. 3.3.3) is con- 
firmed by the inspection of the Wavelet Planes. Four of the 
sources present in WP2 are also detected in WP3 (marked 
CIA, C1B, C1C and C1D in Fig. 3b). The first two are 
still detected in WP4, where the smoothing is so strong 
that only an extended region remains, comparable to CI, 
with two local maxima (at CIA and C1B) and an evi- 
dence of distortion by C1C. Note that sources CIA and 
C1B were detected at ~ 3.5 a and cut using the local de- 
tection algorithm DETECT provided in the PSPC Extended 
Source Cookbook. Thus, the excess photons measured in 
Sect. 3.3.2 do not take these sources into account. But, it 
is very likely that the majority of the diffuse gas is related 
to these sources. 

We now examine sources CIA, C1B, and C1C and 
compare with their optical and radio counterparts: 

Source CIA: 

Fig. 4a shows a zoomed optical image of CIA, to- 
gether with WP3 contours (solid) and radio con- 
tours (dashed). There is an NVSS radio source less 
than 10" away from this X-ray source. Both sources 
are coincident with two optical sources. Moreover, in 
Fig. 4a, we verify that the western part of the radio 
source is aligned with the elongated X-ray source and 
with the direction defined by the two optical sources. 
The second of these optical sources is an emission 
line galaxy, HCG 16-9, observed spectroscopically by 
RdC3Z. Its recession velocity was measured as more 
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Fig. 4. Zoomed optical images of CIA (a) and C1B (b), to- 
gether with WP3 X-ray (solid) contours and NVSS radio 
(dashed) contours. Note the superposition of the X-ray and ra- 
dio sources and in (a) the quasi-perfect alignment of their axis 
angle with the direction defined by the two optical sources. 
The brighter one, HCG 16-9, is a background galaxy (Ribeiro 
et al. 1996). The scale (upper left) corresponds to the distance 
of the HCG 16 group. 



than 20 000 km s -1 , i.e. it is a background galaxy su- 
perimposed on the group. The other optical source is 
comparable in size and magnitude to HCG 16-9, and 
although we don't have any spectroscopic evidence for 
this, it is unlikely that this galaxy lies at the redshift 
of HCG 16. 

The coincidence of an X-ray, a radio and two op- 
tical sources, with one having a radial velocity 
16 000 kms -1 larger than the group makes it very un- 
likely that CIA is linked to the group. The radio source 
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Fig. 3. Contour map of the wavelet planes WP2 (a), WP3 (b), WP4 (c), and WP5 (d) superimposed on an optical DSS image. 
Contours are drawn at 1 a, 2 a, 3 a, 5 a, 10 a and 50 a, where a 2 is the 3 — a clipped variance of the whole WP (see Sect. 3.4.1). 
Structures in WP2, WP3 and WP4 are resolved at 1, 2, 4, and 8' (23, 46, 92, and 184 kpc), respectively. In Fig. 3b, the 
objects XR1 and XR2, detected by Saracco & Ciliegi (1995) are also shown, as well as the different components of region CI. 



denotes a point source, but the X-ray emission seems to 
be extended. It is difficult to say much more about this 
source. If it is at the distance of HCG 16-9 (z - 0.072), 
its extension is several hundreds of kpc. It may be a 
background X-ray group with two prominent galaxies. 
Source C1B: 

Fig. 4b shows that an NVSS radio source is overlap- 
ping the X-ray extended structure, at 30" from the 
X-ray local maximum. The probability that this su- 



perposition is random is less than 0.01, without taking 
into account the X-ray position uncertainty (see Con- 
don et al. 1998). No optical source corresponding to 
C1B is detected cither in the DSS or in dc Carvalho 
et al. (1994). Here again, the superposition of a radio 
and an X-ray source (with no optical counterpart this 
time) puts serious doubts on the link between C1B and 
the group and on the diffuse nature of CI. Moreover, 
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only 10 ± 8 net counts originate from C1B (after the 
point source is masked). 
Source C1C: 

A total of 33 ± 13 net counts are emitted from re- 
gion C1C, which is detected as a local maximum in 
WP2 and WP3. No radio counterpart is detected, but 
a star is found 17" from the X-ray peak (Figs. 2b and 
3b), whose magnitude is Bj = 16.66 according to the 
ROE/NRL COSMOS UKST Southern Sky Object Cat- 
alog. Given that according to the COSMOS catalog 
there are N = 17 brighter stars in a 16' x 16' box 
centered on HCG 16, Poisson statistics yield a proba- 
bility of P = 1 - exp [itN{6/W) 2 ] = 1.7% of having 1 
or more stars within 17" from the center of C1C. We 
thus conclude that the X-ray emission of C1C is linked 
to a foreground star. The Lx / L opt ratio (as defined by 
Motch et al. 1998) is 0.01, i.e. in the upper range for 
stars (see Fig. 3 of Motch et al.). Note that its X-ray 
emission was not excised because it was just below the 
3 a excision threshold. 
Source C1D: 

15 ± 8 net counts (8% of the counts in regions C2- 
C3-C4, and 16% of the counts in CI, see Table 5) are 
detected in region C1D, which has no radio nor optical 
counterpart. We do not incorporate this source in our 
spectral analysis, but do include it in our total group 
diffuse luminosity. 

In summary, of the four components of CI, one (CIA) 
is linked to a background galaxy and a radio-source, one 
(C1B) is a radio-source, one (C1C) is a foreground star, 
and one (C1D) contributes marginally at best to the group 
diffuse emission. Therefore, at best, only one-sixth of the 
X-ray emission from CI is connected to diffuse gas in 
HCG 16. 

3.4.3. Nature of regions C4 and C5 

The same analysis can be repeated with regions C4 and C5 
(regions C2 and C3 are directly around the galaxies of the 
group, and the excess photons are likely to be produced 
in the group). 

Region C4: 

Region C4 has an extended structure with three promi- 
nent local maxima, even in Fig. 1, where the FWHM 
of the smoothing is only 45". The three aligned max- 
ima are still apparent at Co = 25 (Fig. 2c) and some- 
what at Co = 50. (Fig. 2b). They are also recovered in 
WP2. A double structure remains in WP3 and is finally 
merged in WP4. No radio or identified optical source 
is found superimposed on C4. This reinforces its dif- 
fuse nature, without ensuring that this excess is linked 
to the group, even if this seems a reasonable assump- 
tion. Indeed, the two interacting galaxies HCG 16a and 
16b lie at ~ 140 h^kpc in projected distance from C4. 
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Fig. 5. Radio (dashed) and WP2 (a), WP3 (b) wavelet plane 
(solid) contours of C5 superimposed on an optical image. Note 
the several maxima comprising the X-ray emission, and the 
numerous unidentified optical sources present in this area. 



This gas could have been ejected by these galaxies (see 
Sect. 5). 
— Region C5: 

Figs. 5a and b show the superposition of optical, radio 
and X-ray wavelet planes WP2 and WP3, respectively. 
Region C5 has a complex X-ray structure, and several 
local maxima can be seen in WP2, which merge into 
a single structure in WP3. An 8.2 mJy NVSS radio 
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source is superposed with the extended structure in 
WP3, and the prominent three X-ray sources are situ- 
ated at distances of 35, 80, and 95" from the peak of 
the radio emission. Given cumulative source counts of 
10 6 and 4 x 10 5 for 1.4 GHz radio fluxes greater than 5 
and 15mJy respectively (Condon et al. 1998), and the 
82% NVSS coverage of the celestial sphere, we infer 
20.7 NVSS sources brighter than 8.2 mJy per square 
degree. Therefore, Poisson statistics yield respective 
probabilities of 0.6%, 3.2% and 4.5% of having each 
of the three X-ray sources as close as they are to the 
radio-source. 

Likewise, at least six optical sources are clearly visi- 
ble within the outer WP3 isophotes of C5, of which 
two are COSMOS galaxies brighter than Bj = 20.1, 
and there is an additional COSMOS galaxy lying just 
10" outside of the outer WP3 isophote (to the East). 
Within the 16' x 16' box centered on the group, there 
are 34 galaxies in COSMOS with Bj < 20.1. Poisson 
statistics then yield a 3.6% probability of having as 
least 3 COSMOS galaxies within 10" of the WP3 outer 
contour (in a 5.5 deg 2 region). 

If these 3 galaxies were at the distance of HCG 16, 
their typical separations would be 20 h§ kpc and they 
would all be less luminous than the SMC. We would 
then have a very compact subgroup of very faint dwarfs 
detached from Hickson's original compact group of 4 
bright galaxies. 

The alternative of a background group or cluster ap- 
pears much more plausible, given that subgroups of 
very faint dwarfs have never been discussed, and that 
there is X-ray emission apparently associated with this 
group or cluster. 

Therefore, the concordance of the radio source with 3 
X-ray sources, coupled with the large galaxy surface 
number density in C5, strongly suggests that C5 is not 
associated with diffuse emission from HCG 16. 



4. Spectral analysis 

The spectra are extracted using the QPSPEC task in the 
IRAF/PROS environment. Again, we restrict our analysis 
to the energy range 0.20 — 2.01 keV unless otherwise noted. 
We first compare the shapes of the spectra (before back- 
ground subtraction) of different regions with that of the 
background. 

4-.1- Spectral shapes 

The shapes of the spectra, Ck , in specific regions are com- 
pared to the shape of the background spectrum, B^, renor- 
malized to the total counts of the spectrum of the region 
of interest. For this, we use likelihoods with Poisson statis- 
tics: 



hi 



EC* 



-1 



J2C k \nB k -J2J2 ln J > 



k j=2 



where k is the spectral channel. We resort to sets of 10 000 
Monte-Carlo trials to derive probabilities that the ln£s of 
simulated background spectra drawn from the observed 
background, measured relative to the normalized back- 
ground spectrum are smaller than the observed ln£ b s 
measured between the spectra of a given region and of 
the normalized background. The simulated spectra are 
drawn from a Poisson distribution to match the normal- 
ized background spectrum in each energy channel, i.e. with 
expected simulated counts B k (£ Ck/J2^k)- 

Table 6. Comparison of spectral shapes with the background 
spectral shape 



Region 


Total counts 


ln£ 


P(ln£ < ln£ obs ) 


C2+C3+C4 


663 


-77.7 


0.026 


C6 


974 


-77.1 


0.258 



The fourth column represents the probability (from Monte- 
Carlo simulations) that the region has a spectral shape consis- 
tent with that of the background. 



3.5. Summary 



We have shown that the detection of diffuse gas within 
a radius of ~ 200 h^kpc around HCG 16 reduces to the 
significant detection of five regions of diffuse X-ray emis- 
sion, filling less than half of the circle. The small-scale 
structure of these regions, together with associations with 
optical and radio point sources, allows us to reject two 
of these regions (CI and C5) as related to point sources 
or background extended sources. There thus remains three 
diffuse emission regions, two (C2 and C3) surrounding the 
four bright galaxies and the other (C4) ~ 140 /i^kpc away 
from the two interacting galaxies HCG 16a and HCG 16b. 
A better understanding of the gas physics in these regions 
requires an examination of their spectra. 



Table 6 shows the resultant probabilities that each 
spectral region has a spectrum with the same shape as 
that of the background. Not only does the region of ex- 
cess emission, C2+C3+C4 present an excess of counts, but 
it also has a spectrum whose shape is significantly differ- 
ent from that of the background, while region C6 with 
negligible net counts has a spectral shape consistent with 
that of the background. Hence, the spectral shape analy- 
sis confirms that the regions of excess counts are indeed 
the locations of X-ray emission and not caused by poor 
background subtraction. 

4-2. Spectral fits 

The background-subtracted, point-source excised spectra 
were then fit to a hydrogen absorbed (Balucinska-Church 
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& McCammon 1992) MEKAL (Mewe et al. 1985, 1986, 
with Fe L calculations by Liedahl et al. 1995 and the ion- 
ization balance from Arnaud & Rothenflug 1985) plasma, 
using XSPEC version 10, with \ 2 minimization. A vi- 
gnetting correction of the background was performed be- 
fore the background-subtracted spectrum was analyzed. 

Table 7 shows the parameters we used for extracting 
the galaxy spectra. Table 8 presents the results of our 
spectral fitting to different regions within HCG 16. There 
remains many instances where energy channels have fewer 
than 10 net counts, and the Poisson statistics do not re- 
semble gaussians, hence our \ 2 spectral fits are not fully 
appropriate. Therefore, the reduced \ 2 values and the 90% 
error bars given in Table 8 should not be over-interpreted. 



Table 7. Extraction parameters for galaxy spectra 



Galaxy 




RA 


Dec 




radius 








(J2000) 




(arcmin) 


HCG 16a&b 


2 l 


1 09 m 24 ! :0 


-10°07 


49" 


1.5 


HCG 16c 


2 1 


'09 m 38!2 


-10°08 


49" 


1.2 


HCG 16d 


2 1 


'09 m 43"9 


-10°10 


58" 


1.1 


HCG 16-3 


2 1 


'10 m 16=5 


-10°19 


11" 


1.1-0.9 



The spectrum of HCG 16-3 was extracted from an ellipse. 



The low net counts in our background-subtracted spec- 
tra (Table 8) make it difficult to perform reliable spectral 
fits for temperature, metal abundance, and absorbing col- 
umn. Indeed, decent spectral fits require at least 500 net 
counts, whereas we have typically 2 to 7 times less. Never- 
theless, the spectra often do have sufficient counts to pro- 
vide decent constraints on the gas temperature, as well as 
on the bolometric luminosity, once metal abundance and 
absorbing Hydrogen column are fixed to reasonable values. 
Note the large absorption corrections on the bolometric 
luminosity, especially for the lower best fit temperatures, 
where the emission is peaked at low energies, which are 
the most seriously affected by absorption. 

The temperature of the diffuse emission is well con- 
strained: the diffuse emission is cool at kT = 0.27t° 28 
(with much narrower 90% confidence intervals for fixed 
metal abundance). All fits produce an upper- limit (90% 
confidence level) of 0.56 keV. The metallicity is completely 
unconstrained. 

The spectral fits from region C4 (offset from the group 
galaxies) are similar to that of the total diffuse emission, 
and C4 accounts for roughly half of the total luminosity 
of the diffuse emission. However, C4 could be hot if its 
metal abundance is low. Fig. 6 shows the best-fit MEKAL 
spectrum for region C2+C3+C4 together with the residual 
error per bin. The fit is adequate, hence no additional 
galactic absorption nor additional component is required 
by the spectrum. 

The galaxy pair HCG 16a&b and galaxy HCG 16c 
both emit more X-ray luminosity, before correction for 
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Fig. 6. Spectral fits to the diffuse emission in HCG 16 (re- 
gions C2+C3+C4). The solid and dashed histograms repre- 
sent the best fit MEKAL models for solar and one-tenth solar 
metallicity plasmas, with a galactic absorbing column (N H — 
2.0 x 10 2() cm" 2 ). 



absorption, than the group, illustrating the difficulty of 
absorption corrections in low temperature spectra. 

With only 83 photons the spectrum of HCG 16d pro- 
duces virtually no constraints on temperature or metal- 
licity. Nevertheless, careful inspection of the spectrum re- 
veals that the 6 first energy channels (between 0.2 and 
0.5 keV) have very low net counts, even taking into ac- 
count the great error bars, suggesting strong hydrogen 
absorption, as confirmed by the best fits with variable 
absorption, which produce a column density 20 times the 
galactic value. However, the absorption correction to lumi- 
nosity then becomes enormous and produces unreasonably 
high absorption-corrected bolometric luminosities. 

4-3. Comparison with other X-ray results on HCG 16 
4.3.1. Group diffuse emission 

The EINSTEIN satellite pointed for 3217s at HCG 16 with 
the IPC detector. The analysis of these observations by 
Bahcall et al. (1984) led to 26 ± 12 counts for the entire 
group (galaxies plus intcrgalactic medium), corresponding 
to a luminosity of 2 x 10 41 h^ ergs -1 . But the limited 
sensitivity and angular resolution of the IPC didn't allow 
separating the galactic and intergalactic components. 

Saracco & Ciliegi's (1995) analysis of the ROSAT 
PSPC observation of HCG 16 led to the conclusion that 
the X-ray emission was due to point sources associated to 
the galaxies and not to a diffuse intra-group medium. By 
assuming kT — 1.0 keV and Z = 0.1 Zq, they obtained a 



5 a upper limit on diffuse emission of 3.0 x 10 40 erg s _1 
in the 0.5 — 2.3 keV band, after converting this flux to lu- 
minosity correctly (see footnote in Sect. 1 and Table 8), 
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Table 8. Spectral fits 



Region 


Energy range 


Net 


kT 


Z 


N H 


h1 VEM 


x 2 


,2 r lO.B-S.3J 


l2 rune 
"50 ^bol 


^50 Lbol 




(keV) 


Counts 


(keV) 


(Ze) 


(10 20 cm- 2 ) 


(10 63 cm- 3 ) 


(/d.o.f.) 


(10 4J 


ergs -1 ) 




C2+C3+C4 


0.20-1.40 


125 


n 07+0.28 
24+ - 17 

u -^-0.07 
1Q+ 011 

0.19±°; o g 


(0.05) 
0.1 
1 

10 


galactic 
galactic 
galactic 
galactic 


31.4 
24.5 
4.5 
0.5 


0.64 
0.61 
0.69 
0.71 


0.39 
0.43 
0.48 
0.48 


0.68 
0.69 
0.71 
0.72 


2.13 
2.32 
3.31 
3.69 


C4 


0.20-1.31 


71 


0.26tj£I 
20+ - 27 

u - zu -0.07 

on+ - 22 


0.1 

1 

10 


galactic 
galactic 
galactic 


11.9 
2.3 
0.3 


0.30 
0.33 
0.33 


0.23 
0.26 
0.27 


0.36 
0.38 
0.38 


1.03 
1.53 
1.69 


HCG 16a&b 


0.20-2.01 


232 


n 72+ - 24 

fl 79+0.25 
u - ' z -0.18 

0.70±° ? 3 


17+ - 43 
0.1 
1 


galactic 
galactic 
galactic 


16.5 
21.8 
4.0 


0.58 
0.63 
0.74 


0.97 
0.96 
0.90 


1.21 
1.27 
0.99 


1.99 
2.30 
1.40 


HCG 16c 


0.20-2.01 


213 


0-53l°;i« 
0.651°;^ 
0.53 +( "I 

^ mtJKJ — 0.14 


f| 7q+°° 

U - ' y -0.55 

0.1 
1 


galactic 
galactic 
galactic 


7.3 
20.3 
4.0 


0.39 
0.90 
0.37 


0.93 
0.86 
0.94 


1.01 
1.13 
1.02 


1.42 
2.00 
1.42 


HCG 16d 


0.20-2.01 


83 


0-67lS° 31 
0.39tS° 14 
(0.12) 
(0.11) 


0.1 

1 

0.1 
1 


galactic 
galactic 
(41) 
(44) 


7.2 
1.5 
(123) 
(240) 


0.32 
0.27 
0.19 
0.19 


0.31 
0.31 
0.14 
0.14 


0.41 
0.35 
0.16 
0.16 


0.71 
0.53 
190 
1770 


HCG 16-3 


0.20-2.01 


32 


(0.49) 
(0.29) 


0.1 
1 


galactic 
galactic 


2.3 
0.6 


0.23 
0.21 


0.08 
0.10 


0.11 
0.12 


0.21 
0.23 



Column (1): region. Column (2): energy range used for the spectral fit. Column (3): net counts in the given energy band (hence 
lower than the net counts in Table 5). Column (4): temperature. Column (5): metal abundance. Column (6): absorbing hydrogen 
column density (the galactic value is Nh = 2.02 x 10 20 cm -2 , Stark et al. 1992). Column (7): volume emission measure, defined 
as f y n e n p dV. Column (8): reduced \ 2 (P er degrees of freedom) of fit. Their values are low because the noise is not gaussian. 
Column (9): X-ray luminosity in the 0.5 — 2.3 keV band. Column (10): bolometric X-ray luminosity, uncorrected for absorption. 
Column (11): bolometric X-ray luminosity, corrected for absorption. The error bars are 90% confidence levels for one interesting 
parameter. Values in parentheses indicate provide the best fit when the fit was unconstrained (i.e., when the 90% confidence 
levels could not be determined). Values of the metal abundance without error bars nor parentheses were frozen in the fit. 



i.e., ~ 30% below our fitted luminosity in the same energy 
band. 3 

On the other hand, PBEB, who analyzed the same 
data as Saracco & Ciliegi, but in the full 0.1 — 2.4 keV 
PSPC energy band, found an excess of photons in an 8' 
radius circle corresponding to a bolometric luminosity of 
4.8 x 10 41 h^Q ergs -1 . Their factor of 16 discrepancy with 
Saracco & Ciliegi can be subdivided into three terms: 1) a 
factor 3.4 caused by the correction for Hydrogen (and He- 
lium) absorption, 2) a factor 1.6 caused by the narrower 
energy band used by Saracco & Ciliegi, and 3) a factor 3.0 
unexplained residual term. The adaptive smoothing of the 
image (PBEB's Fig. 2) indeed showed two diffuse exten- 
sions on both sides of the group galaxies, roughly corre- 
sponding to the regions CI and C4 defined in Sect. 3.3.1. 

Our spatial analysis is performed over an intermedi- 
ate energy band ([0.2, 2.0] keV, see Sect. 2.1). From our 
spectral fit to the diffuse emission (regions C2+C3+C4), 



3 adopting kT = 0.3 keV instead of 1.0 keV produces roughly 
the same upper limit. 



we find L x [0.5 - 2.3kcv] = 3.9 x 10 40 h^ 2 ergs -1 , only 
slightly higher than the upper limit of Saracco & Ciliegi. 
Our bolometric luminosity is Lbol = 2.3 x 10 41 ergs -1 
(after adding 9% extra luminosity to include region C1D, 
assuming that it has the same spectrum as the rest of 
the group emission) . Thus our derived bolometric luminos- 
ity for the diffuse emission is half that of Ponman et al. 
(1996). 

Summing up from our Table 5, our 8' circle around 
HCG 16 has 306 net counts. In comparison, we infer from 
the surface brightness profile of PBEB (their Fig. 1) that 
they measure 281 net counts within a circle of radius 
0?14 = 8.'4 around HCG 16. In both cases, the group 
galaxies and point sources exceeding a 3cr threshold were 
excised. Thus, the difference (although within the error 
bars) between PBEB's net counts and ours can be at- 
tributed to a different background region. Indeed, their 
background region was chosen to be an annulus of radii 
between 36' and 42', while ours was taken between 26' 
and 34' (see Table 1). We showed in Sect. 2.2 that the 
mean background value increases with off-center radius, 
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which explains our higher net counts. But this can cer- 
tainly not explain the discrepancy in luminosity. Hence, 
the discrepancy between the high luminosity measured by 
Ponman et al. (1996) and our low luminosity for the dif- 
fuse emission is not caused by different background sub- 
traction or different removal of the emission from bright 
galaxies, but by Ponman et al. 's inclusion of regions as- 
sociated with radio or optical galaxies, some of which are 
definitely background objects. 

We do agree with PBEB that this diffuse emission is 
cool, as they find kT = 0.3 ± 0.05(1 a) keV, consistent 
with our values. Note that the correction for absorption 
is as large as 3.5 at temperatures as low as 0.3 keV and 
increases rapidly for decreasing temperature. Hence, while 
our bolometric luminosity, uncorrected for absorption, de- 
pends little on the fit, the absorption corrected luminosity 
is less certain. 

From Table 5, region C6 has 10 ±31 net counts. Hence 
the 90% confidence (1.28 a) and 3 a confidence upper lim- 
its for diffuse emission in C6 are 50 and 103 net counts, 
respectively. We convert these net counts to luminosi- 
ties assuming the same spectrum as region C2+C3+C4, 
which has 172 net counts according to Table 5. After cor- 
rection for galactic absorption, this yields L^oi 

< 6.1 x 

lO^/i^ergs- 1 (90%) and 1.3 x lO 41 ^ 2 ergs" 1 (3a). 
Hence, region C6 contributes little to the diffuse X-ray 
emission of the group (unless it is much cooler, hence more 
absorbed). 



4.3.2. X-ray emission from the group galaxies 

Our 0.5 — 2.3 keV luminosities for the galaxies match fairly 
well those of SC95, who had simply assumed the galaxies 
to have IkeV temperature and solar metallicity: our lu- 
minosities are 1%, 12% and 28% lower for HCG 16a&b, 
HCG 16c, and HCG 16d respectively. 

The extent of the diffuse X-ray emission around each 
galaxy or pair in HCG 16 is fairly large (~ 40 h^ kpc) in 
comparison with the extent of the X-ray emission detected 
by Hcnriksen & Cousineau (1999) around spiral pairs or 
by Read & Ponman (1998) around most interacting pairs. 
Fig. 15 of Read & Ponman shows that only the pairs clos- 
est to maximum interaction show such large extents of 
their diffuse X-ray emission, but their diffuse- X-ray to op- 
tical luminosity ratio are factors of 5 larger than for the 
HCG 16ab pair. 



5. Discussion 

We estimate below the total dynamical mass of the group, 
the baryonic fraction together with the separation of the 
baryonic mass into its different components (namely, hot 
gas, HI, H 2 and stars) and the possibility of virialization 
of HCG 16. 



5.1. Mass budget of HCG 16 
5.1.1. Total dynamical mass 

The total mass of the group can be estimated, to first 
order, assuming that the group is in dynamical equilibrium 
(we will return below to the relevance of this assumption). 
Table 9 presents the estimates for the total mass of the 
group using cither the 4 original galaxies (Hickson 1982) of 
HCG 16, or adding to them the three additional galaxies 
that were found by de Carvalho et al. (1994) in the close 
environment of HCG 16, and confirmed spectroscopically 
by RdC3Z. We apply the mass estimates of Heisler et al. 
(1985) relevant to self-gravitating systems as well as the 
projected mass of Bahcall & Tremaine (1981), relevant to 
test objects orbiting within an underlying potential. We 
use the radial velocity Vi measurements and errors Svi 
from RdC3Z. For the velocity dispersion, we computed 



N 



1/2 



(3) 



We adopt the median of the four Heisler et al. mass 
estimates, as we deem it unlikely that the galaxies are 
test particles in a potential, since we do not detect this 
potential in diffuse X-rays. Interpolating the total mass 
at 8', within which we have constraints on the gas content 
of the group, then yields M to t — 1-4 x \Q 12 h^M Q . For 
a spherical group, this yields a density that is 256 times 
the critical density of the Universe or 770 times the mean 
density if £!o = 0.3. 

Of course, the reliable estimation of the total grav- 
itating mass of the group is difficult with only four to 
seven galaxies. Moreover, the group may not be in virial 
equilibrium (see Sect. 5.2.1 below). Mamon (1993, 1995) 
has quantified the effects of departures from virial equilib- 
rium on the estimation of the masses of groups, initially 
following the Hubble expansion, taking into account the 
softened nature of galaxy potentials. If the galaxies were 
point masses, the mass of a galaxy system near full col- 
lapse should be half of the mass inferred from dynami- 
cal equilibrium (hereafter virial mass), as is well known. 
But since galaxies have softened potentials, the velocity at 
closest approach is only a little larger than for the future 
virialized system (before it coalesces). So, the virial mass 
should provide an adequate estimate of the mass — within 
the apparent radius of a galaxy system near full collapse. 
Hence the masses given in Table 9 are probably roughly 
correct, unless projection effects are important in HCG 16. 

5.1.2. Mass in hot gas 

An upper limit to the observed mass of the diffuse inter- 
galactic gaseous medium (IGM) can be estimated from 
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Table 9. Dynamical mass estimates 



Group 


N 


e 


R 


<7v 


Jw vir 


M ■ 


M 




J,/J proj 






(arcmin) 


(Ko k P c ) 


(kms" 1 ) 




{h£ 1O"M 


) 




Hickson (1982) 


4 


3.2 


74 


99 


l.l 


1.1 


0.9 


0.7 


0.4 


Ribeiro et al. (1996) 


7 


16.9 


389 


76 


2.0 


2.0 


2.0 


1.9 


0.8 



Column (2): number of galaxies. Columns (3) and (4): radius of the smallest circumscribed circle containing the N galaxies, 
in arcmin and in kpc respectively. Column (5): unbiased sample velocity dispersion (corrected for the measurement errors, 
see Eq. [3]). Columns (6), (7), (8), and (9): virial, projected, average and median mass estimates from Heisler et al. (1985), 
respectively. Column (10): projected mass estimate for isotropic orbits from Bahcall & Tremaine (1981). These mass estimates 
do not include correction for measurement errors, and are thus slight overestimates. 



the normalization of the MEKAL plasma model, which is 
the Volume Emission Measure (VEM) defined by 



VEM 



n e n p dV 



n e n p V 



(4) 



where V is the volume of the emitting region, n e and n p 
are the electron and proton densities. Including the con- 
tribution from Helium, the gas mass is 



M] 



hot 



m p n p [l + Y/(l-Y)]V 



(5) 



where m p is the proton mass and Y the Helium mass 
fraction. Since n p 2 < nf = nJT^/[l +Y/(2- 2Y)] (the 
equality being reached for a plasma of uniform density), 
Eqs. (4) and (5) lead to 



M hc 



< 



m p V 1 ' 2 V¥M 112 
(l-r) 1 /2(i_y/ 2 )i/2 



= 1.2 x W'h 



-5/2 
50 



1/2 



2 VEM 



' cm 



1/2 



M , (6) 



where we took a distance D = 79.2 Mpc and Y — 0.24, 
and where 9 is the angular radius of the extended emission, 
and 9 g is the angular radius at which the emission from a 
possible group galaxy was cut. 

Eq. (6) is valid for spherical diffuse regions. Whereas, 
C3 and C4 are nearly circular, C2 is close to being com- 
prised of two circles. The VEM of C4 is taken from Table 8, 
while that of C2 and C3 are each taken as half the differ- 
ence of the VEM of C2+C3+C4 in Table 8 minus the VEM 
of C4. Moreover, both circular regions of C2 are assumed 
to have the same VEM. Table 10 presents the estimates 
of the upper limit to the mass in diffuse intergalactic gas 
M^ x for the various regions of diffuse emission. Summing 
up the contributions of the different regions of diffuse X- 
ray emission, we obtain M hot < 2.1 x 10 10 /i^ 5 ^ 2 Mq. 

One does not gain much in attempting to fit a model 
to each diffuse region. Indeed, if the gas density profile is 
n(r) = «o/[l + (9 /9c) 2 ] 1 / 2 , as in the isothermal model 
with = 1/3 (close to the slope found by Ponman & 
Bertram 1993, for HCG62), then Eqs. (4) and (5) lead to 



3^2 x ( x 2 + 1)V2 



sinh 1 x 



^uniform 



c 3/2( 3 



tan 



5)1/2 



Table 10. Upper limits to the observed gas content of diffuse 
regions 



Region 


hlo VEM 





e a 


h b ^ M h -r 


max 
Phot 




(10 63 cm" 3 ) 


(arcmin) 


(10 10 M Q ) 




C2A 


4.9 


1.7 


1.2 


0.24 


212 


C2B 


4.9 


1.7 


1.2 


0.24 


212 


C3 


9.8 


1.9 


1.2 


0.49 


269 


C4 


11.9 


2.0 


0.0 


1.2 


412 


C2+C3+C4 


31.5 






2.2 


309 


C6 (face value) 


2.5 


8.0 


0.0 


4.3 


24 


C6 (90%) 


12.6 


8.0 


0.0 


9.6 


53 


C6 (3<t) 


25.9 


8.0 


0.0 


13.8 


76 



The last column is the upper limit to the density of hot gas in 
units of the critical density of the Universe. 



where x = 9/9 c . M hot /M^ loim is always greater than 
0.87 for 9 C > 0/100. If, on the other hand, one assumes 
0=1 as found by Mulchaey & Zabludoff (1998), Eqs. (4) 
and (5) lead to 



M ho 



/\ //"uniform 
IV2 hot 



24 



1/2. 



x-*' 2 [sinlT 1 x-x/{x 2 + l) 3 / 2 ] 
[tan- 1 x + x/{x 2 +l) - 2x/{x 2 +l) 2 } 1/2 



which falls to 0.47 for 9 C = 9/5 (close to what we infer 
was derived by Mulchaey & Zabludoff, but to 0.25 for 
9 C = 9/10. Hence, the uniform approximation for the mass 
of diffuse gas is probably valid to within a factor of two 
or at worst four. 

An upper limit to the gas mass from region C6 can 
be estimated, assuming that it is a sphere of 200 h§ kpc 
radius, and that its temperature and metal abundance are 
the same as the best fit case for regions 02+03+04. The 
upper limits for 06 are provided in Table 10, with the 
face value counts, the 90% upper limit or the 3 a upper 
limit. The gas mass of 06 is poorly constrained in con- 
trast with its diffuse X-ray luminosity. However, Table 10 
shows that if the undetected region ( C6) has similar tem- 
perature, metallicity and dumpiness as the detected re- 
gions (C2+C3+C4), then its mass-density is at most (3a 
limit) one-quarter that of the detected regions. Similarly, 
its mean X-ray surface brightness is at most one-sixth 
that of C2+C3+C4. Hence, the low emission of C6 is not 
merely a statistical fluctuation in X-ray counts, but in- 



16 



S. Dos Santos & G. A. Mamon: Clumpy diffuse X-ray emission from HCG 16 



dicative of a true underdensity in the distribution of hot 
diffuse gas. 

In summary, we obtain four upper limits to the mass 
in diffuse hot gas within 8' from the group center: Mhot < 

2.1 x 10 10 /i 50 5 ^ 2 (if the 10 counts in C6 are just noise), 

6.4 x 10 w h^^ 2 (taking the 10 counts in C6 at face value), 

1.2 x lO 11 ^ 572 (with 90% upper limit on C6), and 1.6 x 
lO 11 /^ 572 (with 3cr limit on C6). 

5.1.3. Baryonic fraction 

Within an 8' radius from the group center, the molecular 
gas mass is M Ha = 7.4 x 10 10 ^ 2 M Q (Leon et al. 1998), 
while the mass in cold HI gas is slightly less than Mhi = 

4.5 x 10 10 /i^ 2 M Q (Williams 1998, whose VLA map shows 
that a small fraction of the diffuse HI emission extends 
beyond the 4' radius circle around the group). Note that 
the contributions of dust and ionized hydrogen (estimated 
by MPABB), though important in comparison with other 
galaxies, are negligible within the mass budget of HCG 16. 
Thus, if the 10 counts measured in region C6 are just noise 
or taken at face value, the mass in diffuse hot gas is much 
smaller than the mass in cold gas and fairly negligible 
within the total mass budget of HCG 16. 

The fraction of baryons within 8' of the optical cen- 
ter of HCG 16 is 

M* + Mhi + M H2 + M hot 
h M tot 

= (^T £ + - 085 )^o 1 + / ll o t ^o 3/2 , (7) 

(8) 

where M*/Lb is the mean stellar mass-to-light ratio for 
the group and /hot is the mass fraction of the group in 
hot gas, and is 1.5%, 5%, 8% or 11% depending on the 
interpretation of the 10 net counts in region C6 (noise, 
face- value, 90% limit and 3<r limit, respectively). 

Because the galaxies in HCG 16 have probably all un- 
dergone fairly recent bursts of star formation (e.g. Ribeiro 
et al. 1996), their stellar mass to blue luminosity ratios 
are probably much lower than for normal spirals (i.e. 
M*/Lb < 2.5). There are two ways to estimate M*/Lb 
for each of the four bright galaxies in HCG 16. 

First, MPABB mentioned that galaxies HCG 16a and 
16c have rotation velocities consistent with the Tully- 
Fisher (1977, hereafter TF) relation, whereas galaxy 
HCG 16d has a rotational velocity at most half of what 
is expected by the TF relation. This translates to a lu- 
minosity that is at least 16 times larger than expected 
from the TF relation. Assuming that such was also the 
case for galaxy HCG 16b (its low rotation velocity on one 
side is also half of what is expected from the TF relation, 
although its high rotation velocity on the other side is con- 
sistent with the TF relation), this leads to M*/Lb = 1.19 
assuming (M,/iB)„ orma i = 2.5 for the group of 4 galaxies. 



Alternatively, M*/Lb can be constrained by the times 
since the last bursts of star formation in each of the 4 
bright galaxies of HCG 16. These times can be inferred 
from the optical colors of the galaxies or from their X- 
ray properties. Table 11 provides the values of M*/Lb in- 
ferred from the colors of the 4 bright galaxies, using Fioc 
and Rocca-Volmerange's (1997) PEGASE spectral evolu- 
tion model, assuming a single burst, a Rana & Basu (1992) 
initial mass function, and solar metallicity. 



Table 11. Colors, ages and stellar masses of the HCG 16 galax- 
ies 



Galaxy 


B T 


Bf 


B-R {B-R) u 


Age 


M t /L B 


h 50 M* 












(Gyr) 




(1O 1O M ) 


HCG 16a 


12.99 


12.76 


0.82 


0.72 


0.50 


0.39 


3.0 


HCG 16b 


13.74 


13.27 


1.07 


0.86 


0.68 


0.52 


2.5 


HCG 16c 


13.40 


13.10 


1.04 


0.91 


0.71 


0.57 


3.2 


HCG 16d 


13.90 


13.42 


1.29 


1.08 


0.95 


0.75 


3.1 



Column (1): Galaxy name. Columns (2), (3), and (4): Asymp- 
totic blue magnitude before and after correction for Galac- 
tic and internal extinction, and mean color within the fiB — 
24.5mag arcsec -2 isophote, all from Hickson, Kindl, & Auman 
(1989). Column (5): Color corrected for galactic and internal 
reddening using (B - R)° = B - R + 0.54(B§, - B T ). Columns 
(6) and (7): age and M*/Lb from the PEGASE spectro- 
photometric evolution model of Fioc & Rocca-Volmerange 
(1997). Column (8): Mass in stars. 

The stellar masses in Table 11 yield M*/Lb = 0.53 for 
the group of 4 galaxies. With this value of M*/Lb, Eq. (8) 
yields a baryonic fraction of 18% or 21% if the emission 
from C6 is respectively noise or taken at face value, and 
as much as 25% (90% limit on C6) or 28% (3 a limit on 
C6). 

Note that if significant star formation occurred before 
the last burst, the mean colors of the galaxies would be 
redder than with the most recent starburst. Therefore, 
blue colors indicate even more recent starbursts than listed 
in Table 11. This in turn leads to lower M*/Lb and an 
even lower baryonic fraction. Moreover, if M*/Lb < 0.55, 
then there is more mass in cold (HI+H 2 ) gas than in stars. 

5.2. The dynamical state of HCG 16 

The knowledge of the dynamical state of the HCGs is pri- 
mordial to assess the reality of these close associations of 
galaxies on the sky. In paper II, we show that the low ve- 
locity dispersion of HCG 16 is indicative of a non-virialized 
dynamical state, but show that the alternative scenario of 
chance alignments is even less likely for HCG 16. Here, 
we just want to briefly emphasize two points: is the pecu- 
liar X-ray morphology found in HCG 16 compatible with 
virialization, and are the X-ray properties deduced from 
ROSAT data reduced here compatible with chance align- 
ment models? 
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5.2.1. Can a virialized group have an irregular X-ray 
morphology? 

If HCG 16 were in a (nearly) virialized state, the galaxy 
halos should have merged, and thus the global group po- 
tential should be fairly smooth. Moreover, the diffuse gas 
associated with these halos should have also merged, and 
reached hydrostatic equilibrium within this smooth po- 
tential. For example, HCG 62 is the archetype of such a 
virialized group, as its X-ray morphology is smooth and 
the diffuse X-ray emission extends well beyond the HCG 
galaxies (Ponman & Bertram 1993). The significant gas 
density at the group center is attested by the presence 
of a cooling flow, as witnessed by the increase in sur- 
face brightness and the cooler temperature of the inner 
50 /igQ 1 kpc (Ponman & Bertram 1993). 

The X-ray morphology of HCG 16 is very different 
from that of HCG 62. The diffuse X-ray emission is situ- 
ated within ~ 50 h^Q kpc around the galaxies, plus in one 
clump (C4) at 140 kpc from the nearest galaxy. This 
clumpy X-ray morphology strongly suggests that HCG 16 
is far from virialization. 

This argument against virialization supposes that in- 
tergalactic gas, at the time of formation of the group, was 
able to settle in hydrostatic equilibrium in the shallow po- 
tential of HCG 16. This point of view can be challenged if 
the specific entropy of the intergalactic gas at this epoch 
is high enough to prevent the gas from collapsing with the 
dark matter (Ponman, private communication). Indeed, 
assuming the infalling gas has the density of the Universe 
at the epoch of group collapse, there is a maximum tem- 
perature above which the specific entropy of the infalling 
gas will be higher than the specific entropy of the gas set- 
tled in equilibrium in the potential of the group. Since the 
gas entropy cannot decrease (unless the gas radiates), this 
constraint gives a maximum temperature above which the 
gas cannot settle in the dark matter potential. 

This (relatively) high specific entropy intergalactic gas 
still lacks direct observation, but it is indirectly inferred 
from the changes in the X-ray surface brightness profiles 
from clusters to groups (Ponman, Cannon, & Navarro 
1999). Also, the high temperature level is required in the 
interpretation of a number of observational facts. For ex- 
ample, the negative result of the Gunn- Peterson effect im- 
plies that the Universe had an overall reionization (and 
therefore probably reheating) phase before z ~ 5. More- 
over, the steepening of the Lx—T relation from clusters to 
groups (PBEB) has so far only been explained in models 
where the infalling gas was preheated at a temperature 
of ~ 0.5 keV (Cavaliere et al. 1997, 1998; Balogh et al. 
1999). Thus, provided that HCG 16 is forming today, we 
can understand its irregular X-ray morphology as a con- 
sequence of its low virial temperature, which does not ex- 
clude that the underlying dark matter potential is relaxed. 
Finally, hydrodynamical cosmological simulations indicate 



that most of the intergalactic gas is in the 10 5 — 10 7 K tem- 
perature range (Cen & Ostriker 1999). 

If intergalactic gas cannot collapse onto the group, 
then the diffuse gas observed in the group originates pri- 
marily from the galaxies of the group, in the form of shock- 
heated tidally stripped gas or hot winds from collective 
supernovae ejecta. In the latter case, the diffuse hot gas 
should be metal-rich. Our spectral fits do not constrain the 
metallicity of the group, whereas PBEB found Z < 0.17 
for the diffuse gas that they detected. Since half of their 
detected emission arises from the radio sources in regions 
CI and C5, we cannot confirm the low metallicity of the 
HCG 16 diffuse gas. Better signal-to-noise observations 
are required to answer this question, in particular better 
constraints on metallicity and its spatial variation in the 
group. 

Whether or not the group is near virial equilibrium, 
the hot gas appears too clumpy to be itself in hydrostatic 
equilibrium within a nearly spherical potential. 

5.2.2. Chance alignment models 

Difficulties in the understanding of the properties of HCGs 
(in particular their short crossing-times) have led some 
authors to the conclusion that most of these objects are 
chance alignments along the line-of-sight within larger 
structures, namely loose groups (Mamon 1986), clusters 
(Walke & Mamon 1989) or cosmological filaments (Hern- 
quist et al. 1995). 

Loose groups are obviously not dense enough to be 
globally near full collapse (although their cores may have 
already collapsed and formed interacting binaries), and 
are more likely to be in their early phases of collapse (Di- 
aferio et al. 1993; Mamon 1993, 1994, 1995), as is our Local 
Group, or even in the late stages of their Hubble expansion 
before their turnaround (Valtonen & Byrd 1986). Thus, 
one does not expect to see diffuse gas in a compact group 
caused by a chance alignment within a collapsing near- 
spherical loose group. Indeed, if loose groups have diffuse 
intergalactic gas, this gas is too tenuous to be observed. 
Moreover, this gas should remain fairly cold until it re- 
laxes with the group potential, and this re-heating must 
await the virialization of the group. Therefore, the gas that 
one expects to observe in X rays within loose groups will 
be associated with the dense substructures within these 
groups, i.e., galaxies and interacting binaries (which are 
expected in chance alignments, Mamon 1992), and with 
gas ejected by supernovae near interacting galaxies (e.g., 
Read & Ponman 1998; Henriksen & Cousineau 1999). 

The X-ray properties of chance alignments within cos- 
mological filaments depend strongly on the dynamical and 
thermal state of the cosmological filament. Such filaments 
appear clearly in cosmological simulations. For example, 
the hydrodynamical simulations of Cen & Ostriker (1999) 
show that most of the intergalactic gas is not only at tem- 
peratures of 10 5 — 10 7 K, but also within filamentary struc- 
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tures. Cen & Ostriker argue that this gas was shock- heated 
mainly by structure collapse and possibly also by super- 
novae. 

It presents a major challenge to detect gas, either 
within filaments of loose groups, at temperatures well be- 
low 0.3 keV, since the gas in these systems should be too 
tenuous to be an efficient emitter or an efficient absorber. 
Moreover, its typical temperature is too cool to observe 
in X-rays, and its detection is difficult in the EUV be- 
cause of contamination from emission from the galactic 
corona and strong absorption from galactic HI. Hence, one 
does not expect to detect widespread diffuse X-ray emis- 
sion in compact groups caused by chance alignments within 
loose groups or cosmological filaments. Therefore, in these 
chance alignments scenarios, the diffuse X-rays detected in 
HCG 16 are associated with interacting pairs of galaxies, 
perhaps emitted by gas that was stripped by tidal interac- 
tions or ejected by galactic winds generated by supernova 
explosions. 

Our revised X-ray luminosity for HCG 16 brings the 
group closer to the extrapolation of the cluster luminosity- 
temperature relation and to the the group Lx — T rela- 
tion of Mulchaey & Zabludoff (1998). Moreover, if there is 
a universal luminosity-temperature relation spanning the 
range from individual galaxies to binaries to real com- 
pact groups such as HCG 62 to rich clusters, then one 
expects that a chance alignment of N equal luminosity 
systems along the line-of-sight will produce a group that 
will be located a factor N in luminosity above this uni- 
versal luminosity-temperature relation. This is consistent 
with the position of HCG 16 in the Lx-T diagram rel- 
ative to the relation of Mulchaey & Zabludoff, but the 
X-ray luminosity of HCG 16 is 300 times too large for its 
low temperature in comparison with PBEB's Lx — T re- 
lation, and chance alignments cannot explain such a large 
luminosity excess. 

In Paper II, we investigate in more detail the possibil- 
ity that HCG 16 occurs as a chance alignment within a 
looser group or a cosmological filament . 

5.3. Concluding remarks 

Due to their low virial temperature (T ~ IkeV), groups 
of galaxies and their X-ray emission are of extreme impor- 
tance in the study of processes affecting the baryonic dif- 
fuse content of systems of galaxies, such as reheating and 
early energy injection (Cavaliere et al. 1997, 1998; Balogh 
et al. 1999; Ponman et al. 1999). Earlier analyses of the 
ROSAT PSPC observation of Hickson Compact Groups 
(Saracco & Ciliegi 1995; Ponman et al. 1996) were based 
on counting photons in a circle surrounding the optical 
center of the group. The present detailed spatial analysis 
of X-ray emission in HCG 16 highlights the importance of 
a multi-wavelength study, in particular for rejecting point 
sources identified in the optical and radio wavebands. The 
clumpy nature of the hot gas in HCG 16 appears quite dif- 



ferent from the more regular diffuse emission seen in other 
compact groups, and HCG 16 loses its previous character- 
istic of being the sole spiral-only compact group with hot 
gas tracing a fairly regular potential well (see Ponman 
et al. 1996; Mulchaey 1999). Our study raises neverthe- 
less the question that other compact groups previously 
detected in X-rays may be significantly contaminated by 
superimposed X-ray (point or extended) sources. 

Another question raised by this study is the inclusion 
of HCG 16 in the Lx — T relation of groups and clus- 
ters. Indeed, this relation measures the trend between dark 
matter total mass and mean baryonic density in systems 
where the diffuse gas is in equilibrium within the dark mat- 
ter potential. Since the diffuse hot gas of HCG 16 appears 
to lie in several (N) clumps, one would then expect that 
the group luminosity be a factor N above the Lx — T 
relation, for given temperature. If PBEB's Lx — T rela- 
tion is correct, then the factor of 300 excess luminosity for 
its temperature in comparison with the extrapolation of 
the PBEB's trend is too large to be consistent with the 
idea that most of the X-ray emission originates from a few 
clumps of hot gas in equilibrium within an underlying po- 
tential. Moreover, in the frame of theories explaining the 
Lx — T relation in groups (Cavaliere et al. 1997, 1998; 
Balogh et al. 1999), gas cannot settle in hydrostatic equi- 
librium into a system with a virial temperature smaller 
than ~ 0.5 keV (see Sect. 5.2.1). These two arguments 
naturally lead to the conclusion that the diffuse gas ob- 
served in HCG 16 originates mostly from the galaxies, 
either through tidal stripping or in galactic winds driven 
by supernova explosions. 

The launch of new X-ray satellites (Chandra) with its 
high spatial and spectral resolution and XMM with its 
greatly enhanced sensitivity and spectral resolution will 
eventually allow to answer these questions. The precise 
nature of the diffuse emission seen in HCG 16 is still dif- 
ficult to constrain, because of the low signal-to-noise of 
the ROSAT/PSPC observations. We expect to pursue our 
study of HCG 16 using archival ASCA data, as well as 
a Chandra observation, which we have obtained on this 
compact group. 
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